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Abstract 

The aim of this paper is to solve the "gift exchange" problem: you are one of n players, and there 
are n wrapped gifts on display; when your turn comes, you can either choose any of the remaining 
wrapped gifts, or you can "steal" a gift from someone who has already unwrapped it, subject to the 
restriction that no gift can be stolen more than a total of a times. The problem is to determine the 
number of ways that the game can be played out, for given values of a and n. Several recurrences 
and explicit formulas are given for these numbers, although some open questions remain. 
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1. The problem 

The following game is sometimes played at parties. A number a (typically 1 or 2) is fixed in 
advance. Each of the n guests brings a wrapped gift, the gifts are placed on a table (this is the 
"pool" of gifts), and slips of paper containing the numbers 1 to n are distributed randomly among 
the guests. The host calls out the numbers 1 through n in order. 

When the number you have been given is called, you can either choose one of the wrapped (and 
so unknown) gifts remaining in the pool, or you can take (or "steal") a gift that some earlier person 
has unwrapped, subject to the restriction that no gift can be "stolen" more than a total of a times. 

If you choose a gift from the pool, you unwrap it and show it to everyone. If a person's gift is 
stolen from them, they immediately get another turn, and can either take a gift from the pool, or 
can steal someone else's gift, subject always to the limit of a thefts per gift. The game ends when 
someone takes the last (nth) gift. 

The problem is to determine the number of possible ways the game can be played out, for given 
values of a and n. 

For example, if a = 1 and n = 3, with guests A, B, C and gifts numbered 1, 2, 3, there are 42 
different scenarios, as follows. We write XN to indicate that guest X took gift iV - it is always clear 
from the context whether the gift was stolen or taken from the pool. Also, provided we multiply 
the final answer by 6, we can assume that the gifts are taken from the pool in the order 1,2,3. 
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There are then seven possibilities: 



A1,B2 
A1,B2 
A1,B2 
A1,B2 
A1,B2 
A1,B1 
A1,B1 



C3 



CI, A3 

C1,A2,B3 

C2,B3 

C2,B1,A3 

A2,C3 

A2,C2,A3 



(1) 



and so the final answer is 6 • 7 = 42. 

If we continue to ignore the factor of n\ due to the order in which the gifts are selected from the 
pool, the number of scenarios for the case a = 1 and n = 1,2, 3, 4, 5 are 1, 2, 7, 37, 266, respectively. 

We noticed that these five terms matched the beginning of entry A001515 in p3], although 
indexed differently. The nth term of A001515 is defined as y n {l), where y n (x) is a Bessel polynomial 
(|S!) [2]) 03])) an d for n = 0,1, 2, 3, 4 the values are 1, 2, 7, 37, 266, respectively. Although there was 
no mention of gift-swapping in that entry, one of the comments there provided enough of a hint to 
lead us to a complete solution of the general problem. 

Comments on the rules 

(i) If a = 1 then once a gift has been stolen it can never be stolen again. 

(ii) If a = 2, and someone steals your gift, then if you wish you may immediately steal it back 
(provided you got it honestly!), and then it cannot be stolen again. Retrieving a gift in this way, 
although permitted by a strict interpretation of the rules, may be prohibited at real parties. 

(iii) A variation of the game allows the last player to take any gift that has been unwrapped, 
regardless of how many times it has already been stolen, as an alternative to taking the last gift 
from the pool. This case only requires minor modifications of the analysis, and we will not consider 
it here. 

(iv) We also ignore the complications caused by the fact that you brought (and wrapped) one of 
the gifts yourself, and so are presumably unlikely to choose it when your number is called. 

2. Connection with partitions of labeled sets 

Let H a (n) be the number of scenarios with n gifts and a limit of a steals, for a > 0, n > 1. Then 
H a (n) is a multiple of n\, and we write H a (n) = n!G CT (n — 1), where in G a (n — 1) we assume that 
the gifts are taken from the pool in the order 1,2, ... ,n. We write n — 1 rather than n as the 
argument of G a because the nth gift plays a special (and less important) role. This also simplifies 
the statement of Theorem [TJ 

In other words, G a (n) is the number of scenarios when there are n + 1 gifts, with a limit of a 
steals per gift, and the gifts are taken from the pool in the order 1, 2, . . . , n + 1. 

As mentioned above, the sequence of values of G\{n) appeared to coincide with entry A001515 
in [14]. One of the interpretations of that sequence (contributed by Robert A. Proctor on April 18, 
2005) involved partitions of a labeled set into blocks, and this was enough of a hint to lead us to 
our first theorem. 
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We recall that the Stirling number of the second kind, S2(i,j), is the number of partitions of 
the labeled set {1, ... , i} into j blocks ([6], [7]), while for h > 1 the /i-restricted Stirling number of 
the second kind, 1S2 (i, j), is the number of partitions of {1, ■ ■ ■ ,i} into j blocks of size at most h 

Theorem 1. For a > and n > 0, 

(o-+l)n 

G a (n)= Y, 4 a+1 \k,n). (2) 

k=n 

Proof. Equation (|2|) is an assertion about G a (n), so we are now discussing scenarios where there 

are n + 1 gifts. For a = 0, Hq(u + 1) = (n + 1)!, so Go(n) = 1, in agreement with 5^ (n, n) = 1- 

We may assume therefore that cr > 1. Let an "action" refer to a player choosing a gift 7, either 
by taking it from the pool or by stealing it from another player. Since we are now assuming that 
the gifts are taken from the pool in order, 7 determines both the player and whether the action 
was to take a gift from the pool or to steal it from another player. So the scenario is fully specified 
simply by the sequence of 7 values, recording which gift is chosen at each action. For example, the 
scenarios in ([TJ are represented by the sequences 123, 1213, 12123, 1223, 12213, 1123, 11223. Since 
the game ends as soon as the (n + l)st gift is selected, the number of actions is at least n + 1 and 
at most (<7 + l)n + 1. 

The sequence of 7 values is therefore a sequence of integers from {1, . . . ,n + 1} which begins 
with 1, ends with n + 1, where each number i€ {l,...,n} appears at least once and at most a + 1 
times and n + 1 appears just once, and in which the first i can appear only after i — 1 has appeared. 
Conversely, any sequence with these properties determines a unique scenario. 

Let k denote the length of the sequence with the last entry (the unique n + 1) deleted. We 
map this shortened sequence to a partition of [1, . . . , k] into n blocks: the first block records the 
positions of the l's, the second block records the positions of the 2's, . . ., and the nth block records 
the positions of the n's. Continuing the example, for the seven sequences above, the values of k 
and the corresponding partitions are as shown in Table 1. 

Table 1: Values of k and partitions corresponding to the scenarios in (JTJ) . 

k partition 

2 1,2 

3 13,2 

4 13,24 

3 1,23 

4 14, 23 

3 12,3 

4 12,34 

The number of such partitions is precisely s!f +1 ^ (k, n). Since the mapping from sequences to 
partitions is completely reversible, the desired result follows. ■ 

Remark. The sums B(i) := Ylj^ihj) are t ne classical Bell numbers. The sums Ylj ^2 (*>j) 
also have a long history [10], [TTJ. However, the sums X/j'S'i (*>j) mentioned in ([2]) do not seem 
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to have studied before. Note that the limits in ([2]) are the natural limits on the summand k, and 
could be omitted. 

To simplify the notation, and to put the most important variable first, let 

E a (n,k):=si a+1 \k,n), (3) 

for a > 0, n > 0, k > 0. In words, E a (n, k) is the number of partitions of {1, . . . , k} into exactly n 
blocks of sizes in the range [1, . . . , a + 1]. 

For n > 0, E a (n,k) is nonzero only for n < k < (a + l)n. To avoid having to worry about 
negative arguments, we define E a (n, k) to be zero if either n or k is negative. Then 

(a+l)n 

G a (n)= ]T E a (n,k). (4) 

k=n 

Stirling numbers of the second kind satisfy many different recurrences and generating functions 
([U Chap. V]), and to a lesser extent this is also true for E a (n,k). We begin with three general 
properties. 

Theorem 2. (i) Suppose a > 1. Then E a (n, k) = for k < n or k > (<r + l)n, and otherwise, for 
n < k < (a + l)n, 

E a (n,k) = J2[ ^ )E a (n-l,k-l-i). (5) 
(ii) Fora>0,n>0,k>0, 

k\ 

(ai,...,a CT+ i) 

where the sum is over all (a + l)-tuples of nonnegative integers (ai, . . . , a CT +i) satisfying 

ai + a 2 + 03 • • • + a CT+ i = n , 
ai+2a 2 +3a 3 h (cr + ^c^+i = k . (7) 

(Hi) The numbers E a (n,k) have the exponential generating function 
00 (<x+l)n fc 



E E K(n,^| r =exp -U + | + --- + 77 TT )t) 

n=0 k=n 



(8) 



Proof, (i) This is an analog of the "vertical" recurrence for the Stirling numbers ([6j Eq. [3c], p. 209]). 
The idea of the proof is to take a partition of [1, ... ,k], remove the block containing k, and renum- 
ber the remaining parts, (ii) Here is the number of blocks of size i in the partition. This follows 
by standard counting arguments (cf. [61 Th. B, p. 205]). (iii) This is an analog of the "vertical" 
generating function for the Stirling numbers ([H Eq. [2b], p. 206]), and follows directly from (i). ■ 

The recurrence in Theorem [2^i) makes it easy to compute as many values of E a (n,k) as one 
wishes. Tables 3 through 7 show the initial values of Ei(n, k) through E§(n, k), and Table 8 gives 
the initial values of G a {n) for a = through 8. 
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3. The case a = 1 

In the case when a gift can be stolen at most once, from Theorem [2] we have the recurrence 

Etin, k) = Et(n - 1, k - 1) + (k - tfE^n — l,k — 2), 
for n < k < 2n, with E\{n, k) = for k < n and k > 2n; the explicit formula 

k\ 

El{n,k) = (2n-k)\ (k-n)l 2*""' 
for n < k < 2n; and the generating function 



oo 2n 



n=0 fc=n 



El (n,k) x- y - = e 



,x{y+y 2 /^) 



(9) 



(10) 



(11) 



It follows from (gj) that 



Gi(n) 



2n 

E 

k=n 

n 

E 

i=0 



fc! 



(2n - k)\ (k-n)\ 2 k ~ n 
(n + i)l 



(n-i)\ i\ 2 i ' 



(12) 



Equation (|12p shows that the sequence Gi(n) is indeed given by entry A001515 in [14j . That 
entry gives (mostly without proof) several other properties of these numbers, taken from various 
sources, notably Grosswald [8]. We collect some of these properties in the next theorem. Prop- 
erty (hi) is especially interesting, since the following sections will be concerned with attempts to 
generalize it to larger values of a. We recall from [8] that the Bessel polynomial y n (z) is given by 



'M(z) := 



i=0 



(n + i)\z l 
in - i)\ i\ 2 i 



(13) 



Also 2-Po an d (later) iF\ denote hypergeometric functions. 
Theorem 3. (i) 



(ii) 



(Hi) 



G 1 {n) = y n {\) 
G 1 (n) = 2 F 



n + 1, — n i 

_ > _ 2 



Gi(n) = (2n - l)Gi(n -1) + G x {n- 2) 

forn> 2, with Gi(0) = l,Gi(l) = 2. 

(iv) 



n=0 



,l-v / T := 2x 



VI - 2x 



(v) 



r> t \ e ( 2re ) ! 



(14) 
(15) 
(16) 

(17) 
(18) 
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Proof, (i) and (ii) are immediate consequences of f)12[) . 

(iii) We give three proofs of (|16p , (First proof.) Equation (|16p follows from one of the recurrences 
for Bessel polynomials ([SJ Eq. (7), p. 18], [9]). (Second proof.) Alternatively, it is easy to verify 
from {10]) that 

Ei(n,k) = (2n-l)Ei(n-l,k-2)+E 1 (n-2,k-2). (19) 

Our conventions about negative arguments make it unnecessary to put any restrictions on the range 
over which Q19|) holds. By summing (|19j) on A; we obtain (I16p . (Third proof.) The third proof is 
combinatorial. We will show the equivalent statement that for n > 3, 

Gi(n) = Gi(n - 2) + Gi(n - 1) + 2(n - l)Gi(n - 1) . (20) 

We can build a partition counted in G\{n) in three ways. (A) Take a partition P into n — 2 parts 
and adjoin two parts of size 1, {x} and {y}, say, where x, y are elements not in P. This gives 
G\{n — 2) partitions. (B) Take a partition P into n — 1 parts and adjoin a part {x,y} of size 2. 
This gives Ca(n — 1) partitions. (C) Let P be a partition into n — 1 parts and let S 1 be one of the 
parts. If 5 = {u} is a singleton, then 

P\SU{u,x}l) {y} and P \ S U {it, y} U {x} 

are two partitions into n parts. If 5 = {it, v} is a pair, then 

F\5U{u,i}U {v, y} and P \ S U {it, y} U {«, x} 

are two partitions into n parts. So in either case the pair P, S gives rise to two partitions into n 
parts. There are n — 1 choices for S, so in all we obtain 2(n — l)Gi(n — 1) partitions. The argument 
is clearly reversible, and so (I20j) and hence (|16|) follow. 

(iv) Let 

oo n 

<?x(x) := E^W^F 

n=0 

2 '3 4 

kJO tXj iJC 

= l + 2x + 7- + 37- + 266 ¥ + --- . 
By multiplying (|16p by x n /n\ and summing on n from 2 to oo we obtain the differential equation 

g»(x) = 3g[(x) + 2 X g'{(x) + g x { x ) . (21) 

Then the right-hand side of (|17p is the unique solution of (|21|) which satisfies <5i(0) = 1, ^i(O) = 2. 

(v) This follows from (112p . since the terms i = n — 1 and i = n dominate the sum (see also [51 
Eq. (1), p. 124]). . 



4. The case a = 2 



In the case when a gift can be stolen at most once, the problem, as we saw in the previous section, 
turned out to be related to the values of Bessel polynomials, and the principal sequence, Gi(n), 
had been studied before. For a > 2, we appear to be in new territory — for one thing, the sequences 
Giin), Gs(n), . . . were not among the 140,000 existing sequences in [14] . 

These sequences can be computed using Theorem [2j From @, © we have: 

(a+l)n fc! 

Ga(n) = 2 Yl ai !a 2 !...a (T+1 !l!^2!^...( fJ + l)!^+i ' (22) 

k=n (ai,...,a CT+ i) 
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where the inner sum is over all (cr + l)-tuples of nonnegative integers (ai, . . . , cv+i) satisfying ([7]). 
This may be rewritten as a sum of multinomial coefficients: 



G a (n) 



1 



(7+1 (7+1 



(7+1 



EE-E 

U=l«2 = l in=l 



*1 + «2 H 1" in 

h, h, ■ ■■ , in 



(23) 



where i r is the size of the rth part. 

We naturally tried to find analogs of the various parts of Theorem [3] that would hold for a > 2. 
Let us begin with the simplest result, the asymptotic behavior. This is directly analogous to 
Theorem [3]^v) . 



Theorem 4. For fixed a > 1, 



e((<r + l)n)! 
Ga{n) ~ n!(cr + l)! w 



(24) 



Sketch of proof. The two terms corresponding to {k = (a + l)n, a CT +i = n, other aj = 0} and 
{k = (cr + l)n — 1, acr+i = n — 1, a a = 1, other a, = 0} dominate the right-hand side of (|22p . and are 
both equal to ((cr + l)n)!/(n!(cr + l)! n ). Dividing the sum by this quantity gives a converging sum, 
in which a subset of terms approach 1 + 1 + 1/2! + 1/3! + ..., while the others vanish as n — > oo. ■ 

Concerning Theorem [3](i) , we do not know if there is a generalization of Bessel polynomials 
whose value gives ([22]) for a > 2. 

As for Theorem [3]^ ii) , there is a relationship with hypergeometric functions in the case cr = 2. 
From (l6l) we have 



L(fc-n)/2j 

S 2 (n,fc) = J] 

c=max{0,fc— 2n} 

h/2j 

~ ^ (n-r7 + c)!(?7-2c)!c!2''- c 3 c ' 



(2n - k + c)\(k - n - 2c) Id 2 k ^ n ~ c 3 c 
k\ 



(25) 



where rj = k — n (this is the "excess" of k over n) 



Theorem 5. (i) Let r\ = k — n. 

If r\ < n then 



E 2 (n,k) 



If r\ > n then 

E 2 (n,k) 

(ii) 

G 2 {n) = 



(n + rj)\ 
~ r)\(n - rjj\2n 

(77 + n)\ 



1F1 



(2n - rj)\(r] - n)!2 n 3 r '- n 



1F1 



- V /2, -r]/2 + 1/2 _ g 
n — 7] + 1 ' 3 



-n + ?y/2, -n + 77/2 + 1/2 
77 — n + 1 



' 3 



(26) 



(27) 



n-l 



r?=0 
2n 

+ E 



V- (71 + 77)! 
^ r?!(n - rj)\2 r i 2 1 

(n + 77)! 



-77/2,-77/2 + 1/2 . g 
n — 77 + 1 ' 3 



n=n 



(2n - 77)1(77 - n)!2 n 3 ? ?- n 



^1 



-n + 77/2,-71 + 77/2 + 1/2 _ § 

77 — n + 1 ' 3 



(28) 
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Proof, (i) follows from (|25p using the standard rules for converting sums of products of factorials 
to hypergeometric functions (cf. PQ), and (ii) follows from @. ■ 

We can now state the main theorem of this section, which gives analogs of (|19p and (|16|) . 



Theorem 6. (i) 

E 2 {n,k) = (9n 2 -9n + 2)E 2 {n - 1, k - 3)/2 - 5E 2 (n - l,k - l)/2 

+ (9n 2 - 36re + 35)£ 2 (n - 2, Jfc - 4)/2 + 6(n - 1)£ 2 (« - 2, Jfc - 3) - 3£ 2 (t7 - 2, fc - 2)/2 



+ 3(2n - 5)E 2 (n - 3, Jfc - 4) + 5£ 2 (n - 3, Jfc - 3)/2 + 5£ 2 (n - 4, Jfc - 4)/2. 



(29) 



G a (n) = (9n 2 - 9n - 3)G 2 (n - l)/2 

+ (9?7 2 - 24n + 20)G 2 (n - 2)/2 

+ (6n-25/2)G 2 (n-3) + 5G 2 (n-4)/2, 

for n>4, with G 2 (0) = 1, G 2 (l) = 3, G 2 (2) = 31, G 2 (3) = 18252. 



(30) 



Proof, (ii) Eq. (|30p follows by summing (|29[) on fe, just as (|16p followed from (|19j) . 

(i) We give two proofs of (f29j) . The first proof uses ([26]) . ([27]) and Gauss's contiguity relations for 
hypergeometric functions ([2[ §2.1.2], |15[ §14.7]). There are nine E 2 (i,j) terms in (|29p . and each of 
them is given by either (|26p or (|27p . depending on the relationship between i and j. This means that 
six separate cases must be considered, according to whether k > 2n + l, k = 2n, 2n — 1, 2n — 2, 2n — 3 
or fc < 2n— 4. We give the details just for the first case, the other cases being very similar. Assuming 
then that k > 2n + 1, (I26p applies to all nine E 2 (i,j) terms in (I29p . Writing r] = k — n as before, 
and replacing the final argument | in the hypergeometric functions by a new variable z, we must 
show that the expression 



(rj + n)\ 



(rj - n)!(2n - ?/)!2 n 3 r '- n 
9n 2 -9n + 2 (rj + n - 3)\ 

2 (r]-n - l)!(2n - ?/)!2 n - 1 3 r '- n - 1 

5 (T) + n-l)\ 

2 (rj-n + l)!(2n - - 2)!2 n ~ 1 3^ n+1 



+ 



9n 2 - 36n + 35 
2 

- 6(n- 1) 



(rj + n- 4)! 



3 

+ - 



3(2n - 5) 



(rj - n)\{2n - r] - 2)!2 n ~ 2 3 r '- n 
(rj + n — 3)! 
(r/ - n + l)!(2n - 77 - 3)!2"- 2 3 T '- n + 1 
(77 + n- 2)! 
2 (r/ - n + 2)!(2n - 77 - 4)!2 ri " 2 3 r '- n + 2 
(77 + n - 4)! 



(77 - n + 2)!(2n - 77 - 5)!2 n ~ 3 3 r '- n + 2 
5 (77 + n — 3)! 

2 (77 - n + 3)!(2n - 77 - 6)!2 n " 3 3 , ?- n + 3 
5 (77 + 71-4)! 

2 (77 - n + 4)!(2n - 77 - 8)!2 n ~ 4 3 r '- n + 4 



2-^1 
2-^1 
2-^1 
2F1 
2F1 
2F1 
2F1 
2F1 
2F1 



77/2-77,77/2-77+1/2 _ 
77 — 77 + 1 

77/2 - 77,?7/2-77+l/2 _ 

rj — n 

77/2 - 77 + 1, 77/2 - 77 + 3/2 

77 — 77 + 2 

77/2 - 77 + 1, 77/2 - 77 + 3/2 

77 — 77 + 1 

77/2 - 77 + 3/2,77/2 - 77 + 2 

77 — 77 + 2 

77/2-n + 2,77/2 - 77 + 5/2 

77 — 77 + 3 

77/2-77 + 5/2,77/2-77 + 3 

77 — 77 + 3 

77/2-77 + 3,77/2-77 + 7/2 

77 — 77 + 4 

77/2-n + 4,77/2 - 77 + 9/2 

77 — 77 + 5 



(31) 
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vanishes when z 



Using Gauss's contiguity relations, the nine hypergeometric functions in (|31|! 



can all be expressed as linear combinations of just two of them. The computer algebra program 
Maple 11 simplifies^ the above expression to 



(rj + n - 4)1(3,2 



324(77 -n + l)!(2n- 77 



2)!2 n 3^- n z 3 (z- l) 3 
+ 



^1 2-^1 



?2 2 



Fi 



' 77/2 - n + 1,77/2 - n + 3/2 
77 — 77 + 2 

??/2 - 71,77/2 - 77 + 1/2 _ 
77 — n + l 



(32) 



where (j>i and (f> 2 are polynomials in z of degrees 6 and 5 respectively, with coefficients which are 
polynomials in n and 77. Since the exact values of §\ and 4> 2 are not important for the argument, 
we relegate them to Tables 9 and 10 in the Appendix. The above expression clearly vanishes for 



= |, which proves the desired result. 
Second proof. Let 

D 2 (n,k,c) 



(2n - k + c)\(k - 77 - 2c)!c! 2 k - n ~ c 3 c 
denote the first summand in (1251). We look for a recurrence of the form 



(33) 



^2^2^2C(r,s,t)D 2 (77 + r,k + s,c + t) = 0, 



(34) 



r=0 s=0 t=0 



where the coefficients C(r, s, t) depend on 77 but not on k or c, with the property that when summed 
on c it collapses to the appropriately shifted version of (|29p . which is: 

E 2 {n + 4, k + 4) - (9n 2 + 63n + llQ>)E 2 {n + 3, k + l)/2 + 5E 2 (n + 3,k + 3)/2 

- (9?7 2 + 36?7 + 35)E 2 (n + 2, fc)/2 + 6(77 + 3)E 2 (n + 2, A; + 1) + 3E 2 (n + 2, A; + 2)/2 

- 3(277 + 3)£ , 2(77 + l,/c) -5^ 2 (n + l,/c + l)/2- 5E 2 (n, fc)/2 = 0. (35) 



For this we used the method of Sister Mary Celine Fasenmyer, exactly as described in §4.1 of |12j . 
A Maple 11 program found that there is a solution to (|34p in which the coefficients C(n,k,c) involve 
five free parameters, and there is a two-parameter solution which collapses to (|35|) when summed 
on c. The simplest solution (obtained from Maple's solution by setting both free parameters to 
zero) is the following. All the C(r,s,t) are zero except for the following 19 terms: 



C(0,0,1) 


= -8, 


C(2,l,l) 


= -9, 


C(0,0,2) 


= 7, 


C(2,l,2) 


= 3, 


C(0,0,3) 


= -3/2, 


C(2,2,l) 


= 6, 


C(1,0,1) 


= -18, 


C(2,2,2) 


= -6, 


C(l,0,2) 


= 15, 


C(2,2,3) 


= 3/2, 


C(l,0,3) 


= -3, 


C(3,l,0) 


= -9, 


C(l,l,l) 


= -4, 


C(3,3,l) 


= 5, 


C(l,l,2) 


= 3/2, 


(7(3,3,2) 


= "5/2, 


C(2,0,0) 


= -9, 


C(4,4,l) 


= 1, 


C(2,0,l) 


= 9. 







1 We don't actually know how Maple obtains (|32[) . but the result is consistent with the use of Gauss's relations. 
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It is easy to verify that this collapses to ([35]) when summed on c. 



Is there a combinatorial proof for (130p ? We do not know. 

We discovered (|30p by experiment, using Theorem [6] to suggest the leading term. (Note that if 
r(n) denotes the right-hand side of (I24h . then r{n)/r(n — 1) = (9n 2 — 9n + 2)/2.) We also discovered 
a second recurrence, which is independent of ()30[) : 

(n - 2)G 2 (n) = n(9n 2 - Tin + 17)G 2 (n - l)/2 
+ (6ra 2 - 15n + 13/2)G 2 (n - 2) 

+ (5ra-5)G 2 (n-3)/2, (36) 

for n > 3, with G 2 (0) = 1, G 2 (l) = 3, G 2 (2) = 31. In view of (i28l) . this is equivalent to a 
complicated identity involving hyper geometric functions. We did not find a proof, but Doron 
Zeilberger has kindly informed us that he was able to derive (I36p by applying the method of 
"creative telescoping" ([12l Chap. 6], [17], [18]) to (|33|) and using a modified version of his Maple 
program "MultiZeilberger" . 



5. The case a > 3 



For a > 3 we have not found any connections between G a (n) and generalized Bessel polynomials or 
hypergeometric functions, and we do not have proofs for the recurrences that we have discovered. 

However, we do know that recurrences for G a (n) and E a (n,k) always exist. This follows from 
Wilf and Zeilberger's Fundamental Theorem for Multivariate Sums ( |12] Theorem 4.5.1], [16j). 



Theorem 7. (i) For a > 1, there is a number 5 > such that E a (n,k) satisfies a recurrence of 
the form 

8 S 

E E C if ( n ) E °( n ~hk-j)=0 for all n , (37) 

i=0 j=0 

where the coefficients C\ ■' (n) are polynomials in n with coefficients depending on i and j. 
(ii) For a > 1, there is a number 5 > such that G a (n) satisfies a recurrence of the form 

s 5 

^y^cf ) (n)G a (n-i) = for alln, (38) 

(G) 

where the coefficients C\ (n) are polynomials in n with coefficients depending on i. 

Proof, (ii) As usual, Eq. (f38j) follows by summing (l37jl on k. (i) We will use the case a = 3 to 
illustration of the proof, the general case being similar. We know from ([6]) that 

Mn,k)= £ db] Jl 6c2¥ , (39) 

where the sum is over all 4-tuples of nonnegative integers (a, b, c, d) satisfying 

a + b + c + d = n , 
a + 2b + 3c + Ad = k. 
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In other words, 



Mh^ = y « (40) 

2 n Z-j (2n-k + c + 2d)\(k-n-2c-3d)\c!d!2 k - c 3 c + d ' 1 ; 

where now the sum is over all values of c and d for which the summand is defined. This summand is a 
"holonomic proper- hypergeometric term" , in the sense of [16] , and it follows from the Fundamental 
Theorem in that paper that E^(n, k) /2 n and hence E%(n, k) satisfies a recurrence of the desired form. 
Similarly, in the general case, we write the summand in E a (n, k) as a function of n, k, a 3 , ... , a CT +i, 
again obtaining a holonomic proper-hypergeometric term. ■ 

We conjecture, but do not have a proof, that a stronger result holds, namely that recurrences 
always exist in which the leading terms cffi and C^ G) are both 1, as in ((HJ), (fl~9l) . (l29j) . (|30l) . (HTj) . 
(|42p and Table 11. (The recurrence guaranteed by Theorem [7] may well look more like (|36p . with 
a nontrivial coefficient on the leading term.) 

For a = 3, 4 and 5, we have found recurrences for E a (n, k) and G a (n) with leading coefficient 1, 
although we do not have proofs that they are correct. The following are our conjectured recurrences 
for G 3 (n) and G 4 (n): 



G 3 (n) = (32ra 3 /3-16n 2 + 10n/3-49/6)G 3 (ra-l) 

+ (48n 3 - 236n 2 + 1157ra/3 - 650/3)G 3 (?t, - 2) 

+ (80ra 3 - 382n 2 + 641n-511)G 3 (n-3)/3 

+ (64n 3 /3 - 218n 2 + 2696n/3 - 7915/6)G 3 (?i - 4) 

+ (56n 2 - 490n + 6853/6)G 3 (n - 5) 

+ (56n- 1703/6)G 3 (n- 6) 

+ 58G 3 (n-7)/3, (41) 



G 4 (n) 



1250 n 3 + 625 n 2 



300 n 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



543)G 4 (n- l)/24 
473075?i + 180003)G 4 (n 



184000 n 3 + 447500 n 2 

- 2546500 n 3 + 7679675 n 2 - 12016800 n + 8048577)G 4 (n - 3)/864 



2)/72 

016800n + 8048577)C 

1304291160 n + 1325759504)G 4 (n 



(625 n 4 
(27500 n 4 - 
(336875 n 4 

(4833125 n 4 - 77581625 n 3 + 476892700 n 
(1700625 n 4 + 28316750 n 3 - 605973450 n 2 + 3123850885 n - 5033477363) G 4 (n 
(2670000 n 4 - 64380500 n 3 + 704577200 n 2 - 3610058445 n + 6818722190)G 4 (n 
(2002500 n 4 - 51976000 n 3 + 517392050 n 2 - 2252744530 n + 3561765885)G 4 (n 
(9078000 n 3 - 209915400 n 2 + 1640828980 n - 4301927039)G 4 (n - 8)/7776 
(5393400 n 2 - 91413680 n + 390747263)G 4 (n - 9)/2592 
(1593990 n - 14522219)G 4 (n - ll)/972 
310343G 4 (n- ll)/648. 



4) /2592 

5) /7776 

6) /7776 

7) /7776 



(42) 



The recurrence for G§{n) is similar but more complicated, and we do not state it here. The 
recurrence for E%{n,k) is given in the Appendix (see Table 11). We also omit the recurrences for 
for E±{n, k) and E^{n, k), which are even more complicated. 

Inspection of these recurrences for a < 5 has led us to some conjectures about their general 
structure. First, if J denotes the "depth" of the recurrence, as in (i37|) . (i38j) . then the initial values of 
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S for both G a {n) and E a {n, k) appear to be as shown in Table 2, that is, it appears that these both 
recurrences have depth 5 = ("J 1 ) + 1 (sequence A000124 of [H]). Second, we make the following 



Table 2: Depth 5 of recurrences for G a (n) and E a (n,k). 



a 





1 


2 


3 


4 


5 


5 


1 


2 


4 


7 


11 


16 



conjecture^! about the coefficients in the putative recurrence for E a (n, k). We write this recurrence 



as 



EE^?( n )^( n -^-^ = ' ( 43 ) 

i=0 j=0 

where 5 = ( n ^ X ) + 1, (n) = 1. Then we believe that Cjj(n) = if j > ("J 1 ) + 1, or j < i, or 
(i < a and j > ( n ^ X ) + 1 — ((er + 1 — i) 2 — a — i — l)/2). Furthermore, the degree of Ci t j(n) as a 
polynomial in n is < min{<7, j — i}. 



6. Open questions 

We collect here some of the questions that we have mentioned, (i) The case a = 1 corresponds to 
values of Bessel polynomials; is there a notion of generalized Bessl polynomial that could be applied 
for larger values of a? (ii) The case a = 2 can be described using hyper geometric functions; is there 
a notion of generalized hypergeometric function that could be applied for larger values of cr? (iii) 
Is there a combinatorial proof of (|30p ? (iv) Is the conjecture following Theorem [7] concerning the 
existence of recurrences with leading coefficient 1 true? (v) Find proofs that the recurrences (|4ip 
and (|42p are correct, (vi) Establish the conjectures about the general form of the recurrences for 
G a (n) and E a (n,k) that are mentioned at the end of $5] (this includes question (iv) as a special 
case). 
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Appendix 

This Appendix collects various tables and multi-line formulas that would otherwise have dis- 
rupted the flow of the text. 

Notation: E a (n, k) is the number of partitions of {1, . . . , k} into exactly n blocks of sizes in the 
range [1, . . . , a + 1]. Also G a (n) is the number of scenarios when there are n + 1 gifts, with a limit 
of a steals per gift, and the gifts are taken from the pool in the order 1,2,... ,n + 1. 
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Table 3: Values of Ei(n,k). The array itself appears in several versions in [14] : see for example 
A001498, A144299, A144331; the row sums give A001515, the column sums give A000085 (cf. [TO], 

mi)- 



n\k 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 ... 





1 






































... 


1 





1 


1 
































... 


2 








1 


3 


3 


























... 


3 











1 


6 


15 


15 




















... 


4 














1 


10 


45 


105 


105 














... 


5 

















1 


15 


105 


420 


945 


945 








... 


6 




















1 


21 


210 


1260 


4725 


10395 


10395 


... 


7 























1 


28 


378 


3150 


17325 


62370 


135135 . . . 


8 


























1 


36 


630 


6930 


51975 


945945 . . . 



Table 4: Values of E2(n,k). The array itself appears in several versions in [14] : see for example 
A144385, A144399, A144402; the row sums give A144416, the column sums give A001680 (cf. [JO]). 



n\k 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 ... 





1 






































... 


1 





1 


1 


1 





























... 


2 








1 


3 


7 


10 


10 




















... 


3 











1 


6 


25 


75 


175 


280 


280 











... 


4 














1 


10 


65 


315 


1225 


3780 


9100 


15400 


15400 


... 


5 

















1 


15 


140 


980 


5565 


26145 


102025 


323400 


800800 . . . 


6 




















1 


21 


266 


2520 


19425 


125895 


695695 


3273270 . . . 


7 























1 


28 


462 


5670 


56595 


478170 


3488485 . . . 


8 


























1 


36 


750 


11550 


144375 


1531530 . . . 
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Table 5: Values of E^(n,k). The array itself appears in several versions in [13]: see A144643, 
A144644, A144645; the row sums give A144508, the column sums give A001681 (cf. [TO]). 



n\k 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 ... 





1 






































... 


1 





1 


1 


1 


1 


























... 


2 








1 


3 


7 


15 


25 


35 


35 














... 


3 











1 


6 


25 


90 


280 


770 


1855 


3675 


5775 


5775 


... 


4 














1 


10 


65 


350 


1645 


6930 


26425 


90475 


275275 


725725 . . . 


5 

















1 


15 


140 


1050 


6825 


39795 


211750 


1033725 


4629625 . . . 


6 




















1 


21 


266 


2646 


22575 


172095 


1198120 


7702695 . . . 


7 























1 


28 


462 


5880 


63525 


609840 


5335330 . . . 


8 


























1 


36 


750 


11880 


158235 


1861860 . . . 



Table 6: Values of E^(n,k). The array itself appears in several versions in [13]: see A151338, 
A151509, A151511; the row sums give A144509, the column sums give A110038 (cf. [TU]). 



n\k 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 ... 





1 






































... 


1 





1 


1 


1 


1 


1 























... 


2 








1 


3 


7 


15 


31 


56 


91 


126 


126 








... 


3 











1 


6 


25 


90 


301 


938 


2737 


7455 


18711 


41811 


81081 ... 


4 














1 


10 


65 


350 


1701 


7686 


32725 


132055 


505351 


1824823 . . . 


5 

















1 


15 


140 


1050 


6951 


42315 


241780 


1310925 


6782776 . . . 


6 




















1 


21 


266 


2646 


22827 


179025 


1309000 


9054045 . . . 


7 























1 


28 


462 


5880 


63987 


626472 


5677672 . . . 


8 


























1 


36 


750 


11880 


159027 


1897896 . . . 
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Table 7: Values of E^(n,k). The array itself appears in several versions in [13]: see A151359, 
A151511, A151512; the row sums give A149187, the column sums give A148092 (cf. [TO]). 



n\k 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 ... 





1 






































... 


1 





1 


1 


1 


1 


1 


1 




















... 


2 








1 


3 


7 


15 


31 


63 


119 


210 


336 


462 


462 


... 


3 











1 


6 


25 


90 


301 


966 


2989 


8925 


25641 


70455 


183183 . . . 


4 














1 


10 


65 


350 


1701 


7770 


33985 


143605 


588511 


2341339 . . . 


5 

















1 


15 


140 


1050 


6951 


42525 


246400 


1370985 


7383376 . . . 


6 




















1 


21 


266 


2646 


22827 


179487 


1322860 


9294285 . . . 


7 























1 


28 


462 


5880 


63987 


627396 


5713708 . . . 


8 


























1 


36 


750 


11880 


159027 


1899612 . . . 



Table 8: Number G a {n) of scenarios when there are n + 1 gifts, a limit of a steals per gift and the 
gifts are taken from the pool in the order 1,2, ... ,n. (The array itself is entry A144512 in |14j ; 
the rows give A001515, A144416, A144508, A144509, A149187; the columns A048775, A144511, 
A144662, A147984.) 



n 





1 


2 


3 


4 


5 ... 


G (n) 


1 


1 


1 


1 


1 


1 ... 


Gi(n) 


1 


2 


7 


37 


266 


2431 . . . 


G 2 (n) 


1 


3 


31 


842 


45296 


4061871 . . . 


G 3 (n) 


1 


4 


121 


18252 


7958726 


7528988476 . . . 


G 4 (n) 


1 


5 


456 


405408 


1495388159 


15467641899285 . . . 


G 5 (n) 


1 


6 


1709 


9268549 


295887993624 


34155922905682979 . . . 


G 6 (n) 


1 


7 


6427 


216864652 


60790021361170 


79397199549271412737 . . . 


G 7 (n) 


1 


8 


24301 


5165454442 


12845435390707724 


191739533381111401455478 . . . 


G B (n) 


1 


9 


92368 


124762262630 


2774049143394729653 


476872353039366288373555323 . . . 



16 



Table 9: The polynomial (f>i mentioned in the first proof of Theorem 



</>i = (486 + 729 r/n - 2349 n + 2916 n 2 — 162 77 — 729 n 3 - 729 n 2 rj)z 6 
+ (-45 n + 306 77 - 1080 n 2 + 207 n 3 + 180 77 s - 216 rj 2 - 324 rj 2 n 

- 1539 7777 + 1647 n 2 rj - 54)z 5 

+ (3441 n 3 - 4023 n 2 r) + 24621 n - 348 77 s - 16884 n 2 + 1260 ?? 2 77 

- 10650 7? - 1800 rj 2 + 13203 rjn - 9054)/ 

+ (341 n 3 + 948 r/ 2 n + 8614 77 - 3359 n + 261 n 2 7? + 984 rj 2 

- 6081 r/77 - 484 rf + 3270)z 3 

+ (-35572 77 - 36712 7/77 - 4092 n 3 + 13512 t/ 2 + 2572 t/ 3 + 14952 
+ 25892 7/ + 11164 77 2 7/ - 9244t/ 2 77 + 22472 n 2 )z 2 

+ (11200 t/ 2 77 - 20160 t/ 2 - 3200 t/ 3 - 21120 - 24320 n 2 + 457607/77 
+ 4160 77 3 - 38080 7/ + 42560 77 - 12160 n 2 r])z 

+ (7680 + 7680 n 2 + 14080 7/ - 15360 7/77 - 3840 7? 2 n + 1280 t? 3 

- 14080n - 1280 77 3 + 7680 1] 2 + 3840 n 2 r/) 



Table 10: The polynomial <p2 mentioned in the first proof of Theorem [6l 

(f) 2 = (27 t? 2 + 216 ??77 - 189 1] + 324 + 189 ?7 2 - 675 n)z 5 
+ (-9t] 2 - 9n 2 - 504 t?77 + 495t7 + 9 ?7 + 486)/ 
+ (-15 rj 2 + 600 7777 - 1191 77 - 789 n 2 - 3672 + 3399 n)z 3 
+ (-243 ?? 2 + 408 7777 - 555?? - 303 n 2 - 978 + 1155 n)z 2 
+ (560 ?7 2 - 1360 7777 + 3040 77 - 3440 77 + 720 n 2 + 3840)z 
+ (-320 ?? 2 + 640 7777 - 1600 77 - 320 n 2 - 1920 + 1600 77) . 
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Table 11: The conjectured recurrence for E 3 (n,k). Summing both sides on k gives (|4ip . 



E 3 (n, k) = (32?i 3 /3 - 16n 2 + 22n/3 - l)E 3 (n - 1, k - 4) 

- (4n + 3/2)^3 (n- 1, fc - 2) 

- 17£ 3 (n- l,fc- l)/3 

+ (16n 3 - 88n 2 + 159n - 189/2)£ 3 (n - 2, Jfe - 6) 
+ (32n 3 - 176n 2 + 914n/3 - 497/3)£ 3 (n - 2, jfe - 5) 
+ (28n 2 - 66n + 46)£ 3 (n - 2, jfe - 4) 
+ (-12n + 29/2)£ 3 (n - 2, k - 3) 

- 17£ 3 (n -2, Jfe -2) 

+ (-16n 3 /3 + 152n 2 /3 - 479n/3 + 1001/6)£ 3 (n - 3, jfe - 7) 
+ (32n 3 - 262n 2 + 2218n/3 - 4255/6)£ 3 (n - 3, k - 6) 
+ (84n 2 - 382n + 1247/3)£ 3 (n - 3, k - 5) 
+ (16n - 47/3)£ 3 (n - 3, jfe - 4) 
-28£ 3 (n-3,fc-3) 

+ (64n 3 /3 - 302n 2 + 4154n/3 - 1 2427/6) £3 (n - 4, k - 7) 
+ (84n 2 - 562n + 2858/3)£ 3 (n - 4, Jfe - 6) 
+ (76n - 187)£ 3 (n — 4, — 5) 
-4LE&(n-4,jfe-4)/3 

+ (56n 2 - 574n + 4352/3)£ 3 (n - 5, Jfe - 7) 
+ (84n - 651/2)£ 3 (n - 5, Jfe - 6) 
+ 17£ 3 (n-5,A;-5) 

+ (56n - 1877/6)£ 3 (n - 6, jfe - 7) 
+ 29E 3 (n-6,fe-6) 

+ 58E 3 (n-7,k- 7)/3. 
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